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The spin liquid state of the antiferromagnetic Heisenberg model on a triangular lattice is studied 
within the self-consistent Green's function method. It is shown that the spin excitation spectra is 
gapless, and ground-state energy per site is E g /NJ = —0.966, which is in very good agreement with 
the results obtained within the variational Monte Carlo method based on the resonating-valence- 
bond state. Some thermodynamic properties are also discussed. 
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In the last decades, a large amount of work has been devoted to the understanding of the property of the two- 
dimensional (2D) quantum Heisenberg antiferromagnets. This is motivated partly by the search for quantum dis- 
ordered ground-state and their possible relationship to the high temperature superconductivity Q. The numerical 
simulations strongly support the existence of the antiferromagnetic (AF) long-range-order (AFLRO) for the Heisen- 
berg antiferromagnet on a square lattice with the reduced magnetic moment of about 60% of its classical value [|| . 
It has been shown that this magnetic moment is reduced further by considering the spin frustration arising from the 
next-nearest-neighbors coupling In contrast with the Heisenberg antiferromagnet on the square lattice, the Heisen- 
berg antiferromagnet on a triangular lattice is the three-sublattice spin system with geometric frustration. Anderson 
argued that the strong spin frustration on the triangular lattice may generate a novel resonating-valence-bond 
(RVB) spin liquid ground-state without AFLRO. Since then some different types of quantum ground-state have been 
proposed, in particular, Kalmeyer and Laughlin || argued that RVB state for the triangular lattice is similar to the 
fractional quantum Hall state for bosons. Moreover, various numerical techniques and analytical approaches have been 
used to study this system (|J^. The variational energy reported by Lee and Feng B was essentially identical to the 
result obtained by Anderson ^ based on the RVB state. Huse and Elser Q subsequently constructed wave functions 
exhibiting AFLRO which reportedly had a lower variational energy than Anderson's Q . Therefore the nature of the 
ground-state, especially the existence of AFLRO for the AF Heisenberg antiferromagnet on the triangular lattice, is 
still controversial. 

Apart from the numerical techniques, a popular method to study the quantum spin problem is the spin -wave theory, 
where the quantum spin is mapped onto the boson representation in terms of the Holstein-Primakoff transformation. 
However, the native spin-wave theory, an expansion and linear approximation in Holstein-Primakoff bosons, leads to 
violation of the commutation rule of the quantum spin, this is because that the quantum spins obey the Pauli algebra, 
i.e., the spin one- half raising and lowering operators behave as fermions on the same site and as bosons on different 
sites. In this paper, we study the spin liquid state of the AF Heisenberg antiferromagnet on the triangular lattice 
within the Green's function theory under the Kondo-Yamaji |l(| decoupling scheme. With the help of the spectral 
representations of the correlations, the self-consistent equations are obtained to determine the order parameters. 
The spin excitation spectral is gapless and spin liquid ground-state energy per site is E g /NJ = —0.966. Some 
thermodynamic properties are also discussed. 

The standard AF Heisenberg model on the triangular lattice is written as, 

H = J Sj • (1) 

ir} 

where the sum is over all sites i and, for each i, over nearest-neighbors r\. Since the quantum spin operators obey the 
Pauli spin algebra, then it can be discussed in terms of the Tyablikov's two-time Green's function method [^T). We 
define the spin two-time Green's function G as 

G(i-j,t-t') = -i6(t-t') < [S+(t),Sr(t')} >^« S+(t);S-(t') », (2) 

with 9{t) is the step function, and S~ are the raising and lowering operators of Sf, respectively, and < ■ ■ • > is 
an average over the ensemble. The time-Fourier transform of the two-time Green's function satisfies the equation, 

w « A;B »= — < [A,B] > + « [A,H];B », (3) 
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and the correlation functions can be obtained by the spectral representations as, 



< B(t')A(t) >= i ^ ^ « A;B» u+i0+ -<< A;B»^ l0+ e _ Mt _ tr ^ 
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With the help of the Eq. (3), the equation of motion of the spin Green's function (3) is evaluated as, 

uG(i - j, w) = % < Sf > +2JJ2 « St +v S! - StS? +v ; Sr >>UJ . (5) 

n 

Since the three-sublattice ordering for the triangular lattice is still controversial, and the spin liquid hypothesis is not 
convincingly discarded, then in the following discussions, we only study the system in the case < Sf >= 0, i.e., the 
spin liquid state without AFLRO. In this case, the basic equation for the spin Green's function in the one-dimension 
have been discussed in detail by Kondo and Yamaji [|l0|, and 2D by many authors |Q. Following their discussions, 
we can obtain the spin Green's function (5) as, 

ZJCi 7*, - 1 , , 

G{k,u) = - — , (6) 

where Z is the number of the nearest neighbor sites, 7 /- = (1/Z) e lk ' n , and spin excitation spectra, 

uj 2 (k) = ZJ 2 {± + aC 2 - ad(l + Z lk )]{\ - 7fc ), (7) 

with the order parameters C% =< S^Sf~ +v >, C2 = Ylrrn' iv^v') < ^i+n^i+v' >- ^ n ora - er not to violate the sum rule 
of the correlation function < S^ Sf >— 1/2 in the case < Sf >= 0, the important decoupling parameter a has been 
introduced as these discussed by Kondo and Yamaji |10| ], which is regarded as the vertex corrections. With the help of 
the Green's function (6) and the spectra representation of the correlation functions (4), we obtain the self-consistent 
equations to determine the order parameters Ci,C2, and a as, 

c i = - ZJC ^Y,^^j^ coth[i/Mfc)L (8) 

c 2 = -zjd± £(z 7 2 - coth[~/M*0L (9) 

r-^^w * 1 ^ 4 (10) 

We have performed the numerical calculation for the above self-consistent equations. The ground-state energy per 
site is Eg/NJ = —0.966. For comparison with the other results obtained by the numerical simulations, the present 
theoretical result is also shown in Table I. It is obvious that our result of the ground-state energy is in very good 
agreement with those obtained by Lee and Feng S based on the variational d-wave RVB state without AFLRO, but 
higher than the results of Huse and Elser || for a AFLRO state and Nishimori and Nakanishi's S finite lattice exact 
diagonalization. Huse and Elser || have estimated that the ground-state energy per site is E g /NJ < —1.074, then 
our ground-state energy is more than 10% higher this estimates. But the present low energy of the spin liquid state 
and the short range of its spin correlations indicate that longer-range spin interactions, e.g., second-near-neighbor 
coupling, can push the system further into the liquid phase. More interestingly, the excitation spectra of the RVB 
state for the triangular lattice obtained by Kalmeyer and Laughlin || or Lee and Feng || has a large gap, but in the 
present result, the excitation spectra of the spin liquid state is gapless, which is shown in Fig. 1, and the minimum- 
energy excitation is the spinon at k c = (0,0). Our spin liquid state has spin-i bosonic spinons, and this structure 
clearly disagree with the Kalmeyer and Laughlin [|]|l|] state on the triangular lattice, which was argued to possess 



semionic spinons. We have also computed the spin structure factor S(k) — (l/-Y)X)<ij> < SfSj > e tk '( R * 
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the zero temperature, and the result is plotted in Fig. 2, the maxima of S(k) occur at six wave vectors, which is 
consistent with the result obtained within the large-N expansion p4| |. 

Some interesting properties associated the present spin liquid state, such as, the specific heat and susceptibility 
X = {g 2 ^%/kBT)J2 i j < SfSj > with g is the Landa factor and \ib is the Bohr magneton, have been discussed, 
and the results are shown in Fig. 3(a) and Fig. 3(b), respectively. The specific heat display a broad peak near 
ksT ~ 0.5 J, which is in qualitative agreement with results of the AF Heisenberg model on the square lattice p5[. 
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The susceptibility is very weakly dependent on T, and the similar behavior has been obtained on the square lattice 
fL5f . Therefore, our present results seem to indicate that the global behavior of the AF Heisenberg model on the 
triangular lattice, which is the system with geometric frustration, is qualitative consistent with the frustrated AF 
Heisenberg model on the square lattice. 

In summary, we have studied the spin liquid state of the AF Heisenberg model on the triangular lattice within 
the Kondo-Yamaji's decoupling scheme. It is shown that the spin excitation spectra is gapless, and the ground-state 
energy per site is E g /NJ = —0.966. The energy of this spin liquid is very good consistent with the RVB state, but 
higher than the result for the AFLRO state. The behaviors of the specific heat and susceptibility are similar to these 
obtained on the square lattice. 

A natural question is what is the reason why this self-consistent mean-field Green's function theory is so useful to 
treat the quantum spin systems without AFLRO? To our present understanding, there are at least two reasons: (1) 
the commutation rule (Pauli algebra) of the quantum spin is exactly satisfied in the actual calculations. (2) The sum 
rule of the spin Green's function is always satisfied, and the rotational symmetry is not unphysically broken in this 
self-consistent mean-field Green's function theory. For the spin systems without AFLRO, the low lying excitations 
described are essentially spin waves propagating in a short-range-order with a correlation length. Kondo and Yamaji 
JT(| , and Shimahara and Takada [|l2| have employed this self-consistent mean-field Green's function theory to study the 
one-dimensional Heisenberg spin system and 2D AF Heisenberg spin system on the square lattice, respectively, they 
obtained the results which have satisfactory temperature dependence over whole temperature region from qualitative 
view point. 
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FIG. 1. Momentum dependence of the spin excitation spectra uj(k) of the spin liquid state. The minimum- energy excitation 
is the spinon at k c — (0, 0) 

FIG. 2. Zero temperature structure factor S(k) of the spin liquid state. The global maxima are at six wave vectors. 
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FIG. 3. (a) The specific heat and (b) susceptibility as a function of temperature. 



TABLE I. A comparison of the ground-state energy per site for the antiferromagnetic Heisenberg model on the 
two-dimensional triangular lattice 



Author(s) 


EJNJ 


Method 


Kalmeyer and Laughlin 


-0.94(2) 


Variational RVB 


Huse and Elser 


-1.02(2) 


Variational AF 


Nishimori and Nakanishi 


-1.10(2) 


Finite lattice 


Lee and Feng 


-0.968(4) 


Variational RVB 


The present work 


-0.966 


Green's function method 
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